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ON THE LIMITING STRENGTH OF ROCKS UNDER CON- 
DITIONS OF STRESS EXISTING IN THE 
EARTH'S INTERIOR 



LOUIS VESSOT KING, B.A. (Cantab.) 
Lecturer in Physics, McGill University, Montreal 



§ I. INTRODUCTION 

One of the most important problems in geophysics is to obtain 
a knowledge of the limiting strength of the material forming the 
solid core of the earth as well as that of the principal rocks forming 
its crust. According to the modern theory of elasticity the tend- 
ency of a solid to rupture is measured by the maximum difference 
of the greatest and least principal stresses. This criterion first 
due to Tresca and followed by Sir G. H. Darwin 1 has been found by 
J. J. Guest 2 to give the best agreement with observed results from 
experiments on metal tubes subjected to various systems of com- 
bined stress. Coulomb's suggestion that the greatest shear pro- 
duced in a material is a measure of its tendency to rupture leads 
to a result substantially the same as that given by the hypothesis 
of maximum stress-difference. This follows from the theorem 3 
that the shearing stress at any point is greatest on a plane whose 
normal lies in the plane of the axes of the algebraically greatest 
and algebraically least stress and bisects the angle between them, 
its value being half the algebraic difference between these stress- 
intensities. The family of surfaces whose tangent planes satisfy 
the above condition determines the surfaces along which flaws will 
develop or the material commence to rupture. These surfaces are 
are called surfaces of maximum shear. 

The difficulty in the way of the practical application of these 
criteria lies especially in the determination of the numerical values 
of the limiting stress-difference for different materials. For the 

1 Darwin, "On the Stresses Produced in the Interior of the Earth by the Weight 
of Continents and Mountains," Phil. Trans. Roy. Soc, CLXXIII (1882). 

2 Guest, Phil. Mag. (Ser. 5), XL VIII (1900). 

3 Minchin, Statics, 4th ed., 1889, II, 425. 
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purposes of ordinary engineering practice, a cube of the material 
is put into a testing-machine and the stress required to crush the 
cube is measured. The value of this stress per unit area gives the 
value of the limiting stress-difference for the material as long as 
the boundary conditions are similar to those existing during the test. 
Although this condition generally holds in engineering construc- 
tion work and the criterion is valid, this may by no means be the 
case when the boundary conditions are very different from those 
existing in the test, as for instance in the case of rupture in the 
neighborhood of a cavity or in the case of rupture due to stress- 
difference deep down in the earth's crust. It has been customary 
up to the present to employ the crushing test criterion in applica- 
tions to stresses in the earth's interior. 1 A mathematical analysis 
of the experiments of Dr. Adams, 2 described in the accompanying 
paper, throws a great deal of light on the question and can be made 
to give valuable data on the limiting strength of rocks under con- 
ditions of extreme pressure. 

It will appear that the crushing test criterion of limiting stress- 
difference (e.g., 27,000 lbs. per sq. in. for a cubical specimen of 
Westerly granite of two inches side) is much too low when conditions 
of stress approach those existing in the earth's interior. 

For the purposes of geodynamics, however, "we require to know 
what is the limiting stress-difference under which a material takes 
a permanent set or begins to flow rather than the stress-difference 
under which it breaks; for if the materials of the earth were to 
begin to flow, the continents would sink down and the sea-bottoms 
rise up." 3 

A mathematical analysis of tests described in § 4 afford a valu- 
able means of obtaining information on this question. 

§2. CONDITIONS OF TEST 

The tests described by Dr. Adams on the compression of cylin- 
drical rock specimens containing cylindrical cavities and inclosed 

1 Darwin, loc. cit. 

* Adams, "An Experimental Contribution to the Question of the Depth of the 
Zone of Flow," see this journal, p. 97. 

* Darwin, loc. cit. 
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in heavy nickel-steel jackets lend themselves especially well to a 
determination of stress-differences according to the elastic solid 
theory. An indeterminate condition which might vitiate the 
application of analysis to the problem is the state of initial stress 
brought about by shrinking the nickel-steel jacket over the rock 
specimen. Since the ends of the specimen are free during this 
process it is probable that the rock will adapt itself to a state of 
small initial stress by slight displacements in the direction of the 
axis. This will especially be the case if some time is allowed to 
elapse before the test is carried out. Care was taken to make this 
stress as small as possible under the circumstances; see p. 105. 

One of the most serious objections to the application of the 
theory of elasticity to the determination of conditions of rupture 
lies in the fact that in most tests the conditions of small strains 
upon which the theory is based are violated long before rupture 
takes place. Under the conditions just described the strains remain 
small in consequence of the very slight yielding of the nickel-steel 
jacket. It is well known that the elastic constants are only deter- 
minate so long as the strains are small and that when the material 
is stressed beyond a certain limit the rigidity and compressibility 
take quite different values, so that in most cases the values of the 
stresses just before rupture cannot be calculated from the conditions 
of the test. In the present instance, however, by measuring the 
lateral dilatation of the nickel-steel jacket during a test in the 
manner described in §4, it is possible to determine how the elastic 
constants of the rock change in value with stress. In this way it 
is legitimate to make use of the equations of the elastic solid theory, 
provided we employ instantaneous values of the constants for the 
specimen and as long as the stresses in the nickel steel do not exceed 
the limiting stress for that material. 

§ 3. DETERMINATION OF STRESSES IN ROCK CYLINDER UNDER TEST 

If we consider the state of stress in the rock specimen and in 
the nickel-steel jacket as one of plane stress, the problem may be 
solved without an undue degree of complication. The measure- 
ments given in the next section of the lateral dilatation of the 
nickel-steel jacket when the rock specimen is under pressure show 
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a considerable bulge of the equatorial diameter, showing that the 
problem is not strictly one of plane stress. The complete solution 
of the problem would be of the nature of a similar one considered 
by Filon 1 in which the state of stress is determined for a cylinder 
pressed between planes which are in contact with its ends. For 
the purposes of the present paper, however, we shall deal with the 
problem as one of plane stress and neglect surface tractions, over 
the boundary of the cylindrical specimen due to relative displace- 
ments of the rock surface and the steel jacket when under pressure. 
The results will represent to a first approximation the state of stress 
existing in the rock specimen. From the general solution of the 
problem given by Love 2 it is an easy matter to calculate the stress- 
differences at any point of the specimen. According to the usual 
notation we denote by rr, $$, and zz the radial, transverse, and axial 
stresses respectively: for a case of plane stress there are also the 
principal stresses. We denote the differences of the principal 
stresses by the notation 

(i) fr-zz (ii) rf-U ( iu ) U-zz 

From the theorem cited in § i it is seen that each of the stress- 
differences (i), (ii), (hi) is associated with a family of surfaces of 
maximum shear along which the material will crack or flow. 

(i) The surfaces of maximum principal shear due to the stress- 
differences (fr—zz) consist of two systems of cones of semi-vertical 
angle 45 and cutting each other orthogonally. One such system 
may be imagined to consist of a pile of glass funnels fitting into 
one another. 

(ii) The surfaces of shear due to the stress-differences (fr—OO) 
consist of two systems of mutually orthogonal cylindrical surfaces 
whose traces on a plane perpendicular to the axis of the cylinder 
are equiangular spirals cutting the radii at angles of 45 . 

(iii) The stress-differences (60— zz) give rise to two families of 
helicoidal surfaces of pitch 45 and intersecting orthogonally. The 
traces of these surfaces on the surface of the cylinder give rise to 
helices well known as Luder's fines. 

1 L. N. G. Filon, "On the Elastic Equilibrium of Circular Cylinders under Certain 
Practical Systems of Load," Phil. Trans. Roy. Soc, CXCVIII A, 1902, 182. 
' Love, Elasticity, 2d ed., 140. 
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The systems of surfaces of shear corresponding to each of the 
above cases are shown as (i), (ii), and (iii) in Fig. 1. The cracks 
developed by failure in each of the above ways as they would 
appear on the cylindrical surface, over a cross-section perpendicular 
to the axis and over a section through the axis are represented by 
the diagrams (i), (ii), and (iii). 
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The surfaces of shear developed in a cylinder under test will 
depend on which of the stress-differences is a maximum, a condition 
determined by the elastic constants of the rock and also by the 
nature of the nickel-steel boundary. As soon as this maximum 
stress-difference has attained a limiting value characteristic of the 
material, and determined to some extent by the boundary conditions 
the corresponding families of surfaces of shear will make their 
appearance. 

Considering first the stresses in the rock specimen, the radial 
displacement U satisfies the equation 



3 (W , U\ 
Br\Tr + -;) =0 ' 
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of which the complete primitive is 

U=Ar+j, (i) 

A and B being constants to be determined from the specified bound- 
ary conditions. Longitudinal stress superimposed upon the radial 
stresses can be taken into account by taking the displacement w in 
the direction of the axis to bew=ez, e being a constant representing 
a uniform longitudinal extension. 

If we denote by fr the stress-component along the radius r , by 
$0 the component at right angles to r and by zz that along the 
axis, we find from the stress-equations of equilibrium, 






According to Lame's notation /* is the modulus of rigidity and X 
one of the moduli of elasticity such that *= (A+f m) is the modulus 

\ 
of compression. Poisson's ratio, <r, is denoted by «"= , . . 

In terms of A and B of equation (i) we have for the principal 
stresses and the hydrostatic pressure p 

f}= 2 (X.+ ll .)A-^B+\e 

06=2(k+p)A+^B+\e ) (3) 

zz=2*.A+(\+2it.)e 

p= -^(Fr+ee+zz) = -K(A+e) 

If accented letters refer to corresponding quantities in the 
coaxial nickel-steel jacket, we have for the principal stresses a set 
of equations analogous to (1) and (3), 
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U'=A'r'+B'/r' 

ff'= 2 (X'+ h .')A'-^B'+\'e' 

> (4) 
ee' = 2 (k'+^)A'+^B'+X'e' 

z2' = 2X'^'+(A'+2/x')e' 

Let a be the radius of the cylindrical cavity in the rock specimen, 
b the radius of the specimen, and c the radius of the nickel-steel 
jacket. The boundary conditions give six equations to determine 
the constants A, B, e, A', B', e'. 

At r—a, fr=o, giving 

2(\+ f i)A-2nB/a 2 +\e=o (5) 

At r=r'=b, rr=fr', giving 

2(\+n)A -2i*B/b'+\e=2(k'+p')A'-2 l *.'B'/b>+Xe' (6) 

At r=r'=b, U=U', giving 

Ab+B/b=A'b+B'/b (7) 

At r'=c, fr=o, giving 

2(A'+/»')iT-3/»'£7«H-XV-o (8) 

In addition to these we have in the rock specimen, zz=—P, 
where P is the pressure per unit area applied to the end of the test 
cylinder; this condition gives 

2\A+{\+2tie=-P (9) 

Also since the nickel-steel jacket is free from longitudinal stress, 

zz'=o or 2\'A'+(\.'+2 l i')e'=o (10) 

We denote by ft the quantity, 

I— <T ft' j \ i+<r c'/i—o'/c' /* ) 

We then find from the six equations (5) to (10), 
B = P <r 1 

and } (12) 

A+2ft| 1-20- I+/8 J 
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The principal stress-differences in the rock specimen are from (3) 

(i) fr-zz =—2n{e—A+B/r a ) 

-JL-p J i ~ 2(r t i 1 ; * I 

(ii) fr-fe= -wBlr> =^i> j ^ $ I [> (13) 

(iii) 00-zz=-2/i(e-4-5/r a ) 

= or ( 1-20- ft I_ 0» ) 

The radial displacement U at the outer surface of the rock specimen 
is given by 

ft 2^ i+o- i+ft U4j 

and the longitudinal extension per unit length is given by 

a (i+«r)(i-2ff) 
_P _I_ P-t " I-<r , . 

e ~ 2 M I+ff I+/3 US; 

The radial displacement U' at the outer surface of the nickel-steel 
jacket is given by 

U[_P p i-a*/P «_ , . 

c ~>' c* i-Vl? (i-<t)(i+<t')(i+j8) K } 

and the longitudinal extension per unit length is given by 

e=-o-'— (17) 

c 

§4. DISCUSSION OF EXPERIMENTAL RESULTS 

In order to obtain a numerical verification of the distribution 
of stress in the rock specimens, measurements of the lateral bulging 
of the nickel-steel jacket were made by means of a very sensitive 
Ewing extensometer for a series of increasing loads. The results 
were obtained for a specimen of Solenhofen limestone in which an 
axial and transverse cavity had been drilled, and the experiments 
carried out in the Engineering Testing Laboratory of McGill Uni- 
versity. 
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TABLE I 

Table Giving Radial Displacement of Nickel-Steel Jacket Containing a 
Specimen of Solenhofen Limestone under Varying Loads 



Remarks 





U', -6 


U\, -6 


p 


-X.o 


T Xxo 


Pounds per 


Radial Dis- 


Radial Dis- 


Square Inch 


placement 


placement 




at Center 


over End 


S.ioo 


O 


O 


25,500 


16 


O 


51.000 


34 


4 


76,300 


50 


10 


102,000 


76 


18 


127,000 


106 


24 


153.000 


154 


38 


175.000 


210 


46 


200,000 


326 


66 


5,ioo 


124 


38 


1,000 


124 


38 


1,000 


116 


37-5 



Diameter of nickel-steel jacket 2.5 in., c = i.25; 

length of nickel-steel jacket 3 . 25 in. 
Diameter of rock specimen before loading . 5005 in., 

b =.2502 in.; length of rock specimen before 

loading 1.5725 in. 
Load removed, diam. of specimen at ends .'5052 in.; 

diam. of specimen at equator .5030 in.; length 

1. 5437 in. 
Diameter of cylindrical cavities .050 inch. 
Load removed, vertical cavity filled with rock powder, 

horizontal cavity distorted. 
Longitudinal compression of jacket between 2 in. 

centers .00003 hi., per load of 200,000 pounds. 

Last reading taken 45 hours after preceding 

measurement. 



The nickel-steel jacket used contains about 5 per cent of nickel; 
the elastic constants of nickel-steel containing 55 per cent of nickel 
are taken to be m'= io.8Xio* pounds per sq. in., a J = .327, 1 values 
not very different from these ordinarily given for steel. 

The constants for the Solenhofen limestone of which the rock 
specimen is made have not been determined; the rock in question 
resembles in its structure and properties black Belgian marble for 
which the elastic constants are ^=4.33X10* pounds per sq. in., 
o-=.278. a 

In calculating the value of from (n) we may neglect a 2 /b* 
and since If/c 1 — 04 we find £= . 77. Equation (16) then gives 

^_ P(pound per sq. in.) 
c~ 1. 57X10' 



(18) 



Taking P = 200,000 pounds per sq. in. we have U'/c = 1 2 7 X io -6 . 

The observed radial strain over the center of the specimen is con- 
siderably greater than this. We may attribute the difference to a 
diminution of rigidity. Taking as an extreme case a=o and M=o, 

1 Mercadier, Comptes Rendus, 113, 1801. 

* Adams and Coker, Elastic Constants of Rocks, Carnegie Institution, 1906. 
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o"=5 (incompressible fluid in the nickel-steel jacket), we have 
£=o. Then 

(U'\ _P(poundspersq.in.) . 

which gives for P= 200,000 pounds, 

(7i„-" 8x, °"' 



radial strain 



Vm 



/ 




•opoo. 



POUNDS PER SO. IN. 



Fig. 2. — Showing relation between radial strain and load for a nickel-steel jacket 
containing a specimen of Solenhofen limestone under test. 

Curve I shows the observed radial strain over the center of the specimen. 

Curve II shows the radial strain calculated according to the elastic solid theory 
for rigidity A»=4-33Xio 6 pound per sq. in. and Poisson's ratio <r= . 278. 

Curve III shows the calculated relation between strain and stress for a specimen 
of zero rigidity. 
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The results given in Table I are shown graphically in Fig. 2, 
where the radial strains are plotted on a load base. Equation (16) 
shows that if the rock specimen had behaved throughout as a 
perfectly elastic solid, U'/c should be proportional to the load, and 
the curve connecting load and strain should be a straight line. 
Even taking into account the equatorial bulge of t the nickel-steel 
jacket the proportionality of strain to stress should be exact on the 
elastic solid theory. 1 We may interpret the results just obtained to 
a gradual change in the value of the elastic constants which may be 
regarded as continuous functions of the stress-differences; in fact, 
if the departure of Curve I of Fig. 2 from a straight line be 
attributed to a change in the modulus of rigidity, the variation of 
the latter with stress-difference might be deduced from an analysis 
of the curve. The theoretical curve for the radial strain over the 
center of the specimen agrees well with that observed until a load 
of 80,000 pounds per sq. in. is reached. A curve for the displace- 
ment that would be given if the specimen possessed zero rigidity 
and a Poisson's ratio of \ is also shown in the diagram. The depar- 
ture of the observed curve from that obtained on the hypothesis of 
perfect elasticity shows that the rigidity modulus diminishes under 
load and that Poisson's ratio increases toward the limit \. The 
observations also indicate that when the load is removed the speci- 
men has received a permanent set with a tendency to recovery 
with time, showing an indication of elastico-viscous properties. 2 
The elastic moduli are also permanently changed; a curious fact 
explained by this means is that the specimen after having been 
tested shows a smaller diameter across the center than across the 
ends, whereas the radial displacement was greatest at that point. 
Since the displacements were greatest in the central part of the 
specimen it is reasonable to suppose that the rigidity is permanently 
altered to a lesser value than in the ends. The result is that when 
the load is removed the lateral pressure exerted by the nickel-steel 
jacket is able to produce a greater displacement at the center of the 
specimen than at the ends, so that a slight waist is produced in 
the rock cylinder. 

1 Cf. Filon, loc. cit. 

2 G. H. Darwin, "On the Bodily Tides of Viscous and Semi-elastic Spheroids," 
Phil. Trans. Roy. Soc, CLXX, 17. 



130 LOUIS VESSOT KING 

§ 5. ON THE CONDITIONS OF RUPTURE IN CYCLINDRICAL ROCK SPECI- 
MENS UNDER CONDITIONS OF EXTREME PRESSURE 

An examination of equations (13) shows that the stress- 
differences (i) and (iii) take their maximum values at r—a, and 
that (iii) takes its maximum value at r=&. We thus obtain for 
the maximum values of the stress-differences 

(i) fr— zz=P 
(ii) fr— 



1— o- i-f/3 } (20) 

We note that (i) > (ii) provided /3> 2— !i . Also (i) > (iii) in all 

1 — a 

cases except when fi= °o (b=c) in which case (i) = (iii). 

Thus if 

P>Tz£ ® > (») and (0 > M (") 

so that the specimens will develop the system of surfaces of shear 
shown in (i) Fig. 1. 

If /?< ^^ we have (ii) > (i) > (iii) (2 2) 

In this case the system of surfaces of shear described in (ii), Fig. 1, 
will make their appearance. 
In applying the criterion 

it is necessary to employ instantaneous values of the elastic con- 
stants, since from the previous section these are seen to depend 

on the load. If from (n) we write approximately £=— cor- 
responding to c= 00 , and a=o, the criterion (23) may be written 

As long as <r<% the first of these criteria given in (21) is always 
satisfied. We have seen, however, that the rigidity /* tends to 
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diminish with load, while <r tends to the value |. Under these 

_ it? — 1 

conditions it is quite possible that the condition p<~ given 

1— a 

in (22) may be satisfied before the specimen actually begins to rup- 
ture in the neighborhood of the cylindrical cavity. The break- 
down of a rock specimen under these conditions would be charac- 
terized by a tendency to split across a plane through the axis : over 
a cross-section a system of spiral cracks cutting the cylindrical 
cavity at an angle would make their appearance and long splinters 
parallel to the axis would break off from the interior of the cylin- 
drical cavity. Actual cases of fractures of this kind are described 
by Dr. Adams in the preceding paper. 1 We note that for a solid 
cylinder (a=o) the stress-difference (i) = (iii), so that the conditions 
for the appearance of the surfaces of shear (i) and (iii) are identical. 
System (iii) is characterized by the appearance of Luder's lines over 
the surface of the specimen; these in fact are often observed. 

§ 6. ON THE LIMITING STRESS-DIFFERENCE IN THE NEIGHBORHOOD 
OF A SMALL CYLINDRICAL CAVITY 

The preceding analysis may be applied to a determination of the 
limiting stress-difference S in the neighborhood of the small cylin- 
drical cavities in the specimens of Westerly (Rhode Island) granite 
tested by Dr. Adams and described in the preceding paper. This 
granite gave a crushing strength of 27,370 pounds per sq. in. for 
a two-inch cube. If we note that for this granite <r= .219,* the 

condition p>- is satisfied and the greatest stress-difference 

1— <r 

in the specimen is P, so that S=P. Experiments 373 and 357 

described in Dr. Adams' paper (p. 113) enable a limit to be assigned 

to the value of S. In experiment 373 a pressure P of 160,000 

pounds per sq. in. applied during 75 days produced no change in 

the diameter of the cylindrical cavity, indicating that no flow or 

permanent set had taken place. In the case of a similar specimen 

tested under a pressure P of 200,000 pounds for the same length of 

time, it was found that the vertical cavity was partly closed up by 

1 Adams, p. 109 of this journal, Plate I, figs, c and d; Plate II, fig. c. 

2 Adams and Coker, loc. cit. 
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rock powder. We are therefore able to assert that for a state of 
stress at ordinary temperature and lasting for 75 days, S lies between 
the limits of 160,000 and 200,000 pounds per sq. in., i.e., 

i6o,ooo< S< 200,000 pounds per sq. in. (24) 

Even at a temperature of 55o°C. experiment 381 (p. 115) shows 
that S is greater than 96,000 pounds per sq. in. 

We cannot, however, safely assert that flow or rupture would 
not take place for a stress-difference less than 5 if the boundary 
conditions were very different from these already described at the 
surface of a cylindrical cavity. In order to settle this point experi- 
ments of the type described in § 4 would be required. The value 
given in (24) for S is very much greater than the limiting stress- 
difference of 27,000 pounds per sq. in. obtained by means of the 
usual crushing test. Since this granite is typical of the great mass 
of igneous rocks which make up the earth's crust, it follows that 
in geophysical problems which involve a knowledge of limiting 
stress-difference in the earth's crust, a value considerably higher 
than that usually employed must be taken. 

§ 7. ON THE EXISTENCE OF CAVITIES IN MEDIA UNDER STRESS 

(i) State of Stress in the Neighborhood of Cavities in Elastic Media 

under Shear 
It can be shown that in the neighborhood of both a spherical 
and a cylindrical cavity, the shear can be nearly equal to twice 
the shear at a distance. 1 We see therefore that the medium near 
a cavity will collapse if the stress-difference at a distance exceeds 
half the limiting stress-difference of the rock material. If the 
cavity be small this limiting stress-difference may be taken to be 
the value of S given in equation (24). 

(ii) State of Stress in the Neighborhood of a Spherical Cavity in a 
Medium under Pressure 
If p be the hydrostatic pressure in the medium it can easily be 
shown from the solution given by Love 2 for a spherical shell under 

1 Love, Elasticity, 2d ed., 245. See also a paper by Larmor, Phil. Mag. (Ser. 5), 
XXXIII, 1892, p. 77- 

2 Love, op. tit., 139. 
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external pressure that the maximum stress-difference in the neigh- 
borhood of the cavity is | p. This result is independent of the 
elastic constants and of the radius of the cavity. If then the 
pressure in the medium exceeds f the limiting stress-difference of 
the rock, the cavity will break down. The value of S given in 
(24) for the limiting stress-differences may be employed provided 
the radius of the cavity be small. 

(***) On the Stability of a Spherical Cavity in a Medium under 

Pressure 

That the result quoted in (ii) should be independent of the radius 
of the cavity seems contrary to experience. The explanation of 
this anomaly lies in the fact that the question turns on an exami- 
nation of stability rather than on a purely static consideration 
of stress. A problem of a similar kind is that of the collapse of a 
thin-walled cylinder (e.g. a boiler flue) under external pressure. 1 
The result shows that for a cylindrical shell of infinite length collapse 
will not occur as long as the pressure does not exceed 

where /* and a denote the rigidity modulus and Poisson's ratio for 
the material and t/a is the ratio of the thickness to the diameter. 
The vibrations in the neighborhood of a spherical cavity in an 
infinite medium under hydrostatic pressure do not appear to have 
been explicitly worked out; the problem would correspond to the 
external solution of the problem of the vibrations of a solid elastic 
sphere. 2 In so far as the stability depends on the radial vibrations 
of the cavity, a solution could easily be adapted from Poisson's 
solution for the case of a solid elastic sphere. 3 The result of an 
investigation on stability would give a result analogous to (25) for 
a cylindrical shell. When the pressure p is specified the condition 
of stability will lead to a limiting radius which will be determined 
by the elastic constants of the medium. A cavity having a greater 

1 Love, op. cit., 530. 

' Ci. Love, " Gravitational Stability of the Earth," Phil. Trans., 207 A, 1008. 

> Poisson, "MemoiresurPequilibre et le mouvement des corps elastiques," Mint. 
Paris Acad., t. 8, 1829. 
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radius than this limiting radius will be unstable and collapse in 
such a way as to form a number of smaller cavities each of which 
might be stable. 

§ 8. NOTE ON CONDITIONS OF STRESS IN THE EARTH'S CRUST 

The determination of conditions of stress existing in the earth's 
crust constitutes one of the most difficult problems of geodynamics 
and one on which many recent investigations have been published. 
If the earth is considered as a perfect sphere in which the density 
is a function of the radius only, the problem can be solved without 
difficulty in terms of the elastic solid theory. 1 The application 
of these results throughout the interior of the earth is not legitimate 
because an infinite number of laws of density can be formulated 
to represent the known distribution of mass throughout the earth. 
Near the surface considerations of surface inequalities vitiate the 
result. In this connection Love states, 2 "Apart from the difficulty 
concerning the initial stress in a gravitating body of the size of the 
earth, a difficulty which we seem unable to avoid without treating 
the material as incompressible, there is another difficulty in the 
application of such an analysis to problems concerning compressible 
gravitating bodies. In the analysis we take account of the attrac- 
tion of the inequality at the surface, but we neglect the inequalities 
of the internal attraction which arise from the changes of density 
in the interior; yet these inequalities of attraction are of the same 
order of magnitude as the attraction of the surface inequality." 

In a recent work Love 3 has attacked the problem of isostatic 
support of continents and mountains. "The problem admits of 
an infinite number of solutions even if the distribution of the mass 

is known The problem must remain indeterminate, and 

all we can do is to obtain explicitly one or more of the infinitely 

1 Lam6 solves the problem for an isotropic elastic sphere under its own gravita- 
tion (Love, Elasticity, 2d ed., 140); the result can only be reconciled with fact by the 
assumption of incompressibility or by taking into account a state of initial stress. 
L. M. Hoskins in a paper "Flow and Fracture of Rocks as Related to Structure" 
(U.S. Geological Survey, 1894-95, Part I, 854), solves the problem throughout a com- 
pressible concentric crust of uniform density and small thickness; no account is taken, 
however, of surface inequalities or of considerations of isostasy. 

3 Love, Elasticity, 2d ed., 253. 

3 Love, Some Problems of Geodynamics (Adams Prize Essay, 1911), chaps, ii and 
iii, p. 12. 
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numerous solutions. One solution of the problem was obtained 
by Sir G. Darwin 1 by assuming that the stress is connected with 
a displacement by the same equations as hold in the ordinary theory 
of an elastic incompressible solid, . . . . " 

The theoretical developments turn on the calculation of stress- 
differences corresponding to surface inequalities representing con- 
tinental blocks and mountains. Expressions for the hydrostatic 
pressure at any depth are not explicitly worked out although it 
would be possible to calculate from the analysis given the pressure 
at any depth due to those harmonic inequalities which Love has 
shown can be made to represent the existing distribution of land 
and water and the general form of continental blocks. For 
most purposes it is sufficient to calculate the pressure at any depth 
as that due to the weight of the corresponding column of rock: 
this formula gives a value of greater relative accuracy the greater 
the depth at which the pressure is calculated. 

§9. DEPTH AT WHICH CAVITIES CAN EXIST 

(i) Collapse Due to Stress-Differences in the Earth's Crust 

Love 2 shows that on the isostatic theory the maximum stress- 
difference due to a harmonic inequality of the third order repre- 
senting a continental block of maximum elevation 2 kms. is approxi- 
mately the weight of a column of rock equal to 0.0208 of this 
greatest height, i.e., about 0.0114 of a metric tonne per sq. cm. 
(162 pounds per sq. in.). This value is greater than that brought 
about by inequalities of the first and second order. The stress- 
difference is also worked out for a parallel series of mountain ranges 
400 kilometers apart and having a height of 4 kms. from crest to 
valley-bottom; it is shown that the greatest value of the stress- 
difference is about . 26 metric tonnes per sq. cm. (3,700 pounds per 
sq. in.). In both cases the stress-difference is very much less than 
I S so that no state of stress-difference due to weights of continents 
or mountains is intense enough to cause the material of the 
earth's crust to rupture in the neighborhood of a small cavity. 

1 Darwin, "On the Stresses Caused in the Interior of the Earth by the Weight of 
Continents and Mountains," Phil. Trans. Roy. Soc, CLXXIII (1882). 

2 Love, op. oil., 36; also 47. 
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No account is here taken of stress-difference brought about by 
the folding of the earth's crust or that due to residual effects of 
a state of very intense stress which may have been brought into 
existence at a very remote period of the earth's history. 

(it) Collapse Due to Pressure in the Earth's Crust 
We have noticed in § 8 that the pressure at any point of the 
earth's crust may be roughly represented by the weight of a column 
of rock of that depth. If w represent the average weight of rock 
at the earth's surface per cubic cm. and h the depth in cms., a small 
spherical cavity can exist provided 

wh<%S (26) 

where 5 is the limiting stress-difference of § 6 expressed in grammes 
per sq. cm. Making use of the numerical values corresponding to 
a limiting stress-difference between 160,000 and 200,000 pounds 
per sq. in. (1.13X10 6 and 1.41X10 6 grammes per sq. cm.) for 
Westerly granite and taking 2.7 as the mean density of sur- 
face rocks, we can assert that small cavities in the earth's crust will 
not collapse provided the depth does not exceed a quantity h whose 
value lies between 27.7 and 33.7 kilometers (h between 17.2 and 
20.9 miles). This estimate is several times greater than that 
deduced by Hoskins, 1 the reason being that it is necessary to make 
use of a very much higher value of limiting stress-difference than 
that employed in the investigation referred to. Even at a tem- 
perature of 550° C. which is supposed to exist at a point n miles 
below the earth's surface, Adams has shown that S is greater than 
95,000 pounds per sq. in. which corresponds to a depth of 15 miles. 
Any internal pressure in the cavity (due to water-pressure, steam, 
etc.) would increase the above estimate for the depths at which 
cavities can exist. 

The size of a cavity which can exist at a given depth depends 
on considerations of stability and would demand a separate inves- 
tigation. 

1 Hoskins, op. cit., 850. Hoskins expresses the result of his investigation in the 
form " cavities must certainly close up if wh>S." The value for 5 is taken to be 25,000 
pounds per sq. in. The much higher value between 160,000 and 200,000 pounds per 
sq. in. indicated by experiments on Westerly granite shows that Hoskins' estimate for 
the depth at which cavities must close, must be increased between 7 and 8 times. 
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§ IO. SUMMARY AND CONCLUSIONS 

i. Sees, i to 3 are devoted to a mathematical discussion of the 
state of stress existing in cylindrical rock specimens incased in 
heavy nickel-steel jackets and described by Dr. Adams in an accom- 
panying paper. 

2. In sec. 4 it is shown that observations on the radial strains 
of the nickel-steel jacket agree well with the calculated values as 
long as the load does not exceed a certain value; beyond that 
point the results of experiment indicate a change in the elastic 
constants with load; an analysis of observations of this type could 
be made to give a relation between the rigidity and the load for 
stress-differences beyond the elastic limit. The results of sections 
4 and 5 also give a good account of the manner in which the rock 
specimens break down under load. 

3. The experimental results of Dr. Adams discussed in sec. 6 of 
the present paper show that in the neighborhood of a small cylin- 
drical cavity, in a specimen of Westerly granite, no flow or set 
takes place when the stress-difference amounts to as much as 160,- 
000 pounds per sq. in. The cavity just begins to break down for a 
stress-difference of 200,000 pounds per sq. in. We therefore assign 
to the limiting stress-difference in the neighborhood of small cavities 
in rocks typical of that forming the greater part of the earth's crust, 
a value S between 160,000 and 200,000 pounds per sq. in. for 
ordinary temperatures. 

4. Sees. 7 to 9 deal with the existence of cavities in media under 
stress with especial reference to the depth at which cavities in the 
earth's crust can remain open. 

It is shown that no state of shearing stress in the crust of the 
earth, due to the weights of continents and mountains can cause 
the collapse of the rock in the neighborhood of a small cavity. 

It is also shown that as far as hydrostatic pressure in the earth's 
crust is concerned a small cavity at ordinary temperatures will 
remain open provided the depth does not exceed a value between 
17 . 2 and 20. 9 miles. At a temperature of 550 C. supposed to exist 
1 1 miles below the earth's surface, cavities will remain open when 
submitted to considerably greater pressures than are found at this 
depth. These values greatly exceed previous estimates because 
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experiment shows that a much higher value of limiting stress- 
difference than that usually employed must be taken in the neigh- 
borhood of small cavities. 

5. The size of a cavity which can exist at a given depth depends 
on considerations of stability and would demand a separate investi- 
gation. 

January 30, 191 2 



